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Abstract 

In this paper we investigate the possibility of measuring the general 
relativistic gravitoelectric correction p( GE ) to the orbital period P of 
the transiting exoplanet HD 209458b 'Osiris'. It turns out that the 
predicted magnitude of such an effect is ~ 0.1 s, while the most recent 
determinations of the orbital period of HD 209458b with the photo- 
metric transit method are accurate to ~ 0.01 s. It turns out that the 
major limiting factor is the error in the measurement of the projected 
semiamplitude of the star's radial velocity Km- Indeed, it affects the 
determination of the mass m of the planet which, in turn, induces a 
systematic error in the Keplerian period P^ of ~ 8 s. However, the 
present situation may change if future one-two orders of magnitude 
improvements in the Doppler spectroscopy technique will occur. 

1 Introduction 

In this paper we investigate the possibility of performing a test of general 
relativity, in its weak-field and slow-motion approximation, for the first time 
in a planetary system other than our Solar System 1 . 

More precisely, it has been shown (Soffel 1989; Mashhoon et al. 2001) 
that the general relativistic gravitoelectric (GE) part of the gravitational 
field of a static, spherically symmetric mass 2 also induces an additional small 
correction p( GE ) to the Keplerian orbital period of a test particle. For 

See e.g. http://www.obspm.fr/encycl/ for regularly updated information about the 
extrasolar planets. 

2 It is the Schwarzschild component of the space-time metric tensor which gives rise 
to various tested phenomena (Will 1993) like the well known secular Einstein pericenter 
advance, the deflection of the electromagnetic waves, the Shapiro time delay and the de 
Sitter geodetic precession of the spin of an orbiting gyroscope. 



a circular orbit of radius a around a body of mass M it is, to order O(o 2 ) 

p(GE) „ ^l^GM^, (1) 

where G is the Newtonian gravitational constant and c is the speed of light in 
vacuum. Such an effect has not yet been directly measured. The possibility 
of detecting it in the Sun-Mercury system, for which it amounts to p( GE ) ~ 
0.3 s, has preliminarily been studied in (Iorio 2005). 

At first sight rather surprisingly, it might be easier to detect p( GE ) in 
some of the recently discovered exoplanets than in our Solar System. Indeed, 
many of the exoplanets whose parameters and ephemerides are best known 
have masses comparable to that of Jupiter and orbital periods of just a 
few days. This allows for a high accuracy-to a subsecond level in some 
cases-in measuring their periods from combined photometric transit and 
spectroscopic radial velocity measurements. Indeed, it is more likely to 
observe a transit if the orbital period is short. Moreover, the magnitude of 
the disturbing effects can only be observed if the main star is not too massive 
and the planet is not too light. Of course, such kind of measurements could 
not be performed in the case of Mercury due to its extremely tiny perturbing 
effect on the Sun and its relatively long period of almost 88 days. Only 
difficult and rather sparse astrometric measurements could be done. 

The possibility of testing other relativistic signals which affect the orbital 
motion, like the gravitoelectric secular precessions of the periastron uj and 
of the mean anomaly A4, and the gravitomagentic secular precessions of uj 
and of the longitude of the ascending node £1, must be ruled out because the 
ephemerides of the extrasolar planets do not currently allow to determine 
all the details of the orbital paths except for the orbital period. 

2 The Two-Body Relativistic Correction to the 
Orbital Period 

In the first post-Newtonian approximation, the motion of two particles of 
rest masses M and m, each one moving in the gravitational field of the 
other, can be studied in the post-Newtonian center of mass frame (Soffel 
1989). In the standard isotropic coordinate 3 r, not to be confused with 
the Schwarzschild radial coordinate r which, for M » m, is r = r(l + 

3 It is the coordinate used in the force models of the post-Newtonian equations of motion 
(Estabrook 1971) adopted, e.g., for the computation of the planetary ephemerides by the 
Jet Propulsion Laboratory (JPL). 
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GM/2c 2 r) 2 , the general relativistic equation of motion is (Soffel 1989) 
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^(4 + 2^)-(l + 3^)r 2 + ^(r.r) 2 
r 2H 



r + (r • r)(4 - 2v)r 



(2) 

where [i = G(M + m), v = Mm/(M + m) 2 and the overdot represents the 
time derivative. 

In the case of circular orbits, i.e. r • r = 0, which often occurs for many 
exoplanets orbiting close to their star due to tidal circularization, the post- 
Newtonian acceleration is entirely radial. In polar coordinates @ reduces 
to 
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(3) 



The orbital period can easily be calculated from (J3J) for r =const= a. From 

(4) 
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one gets that the time P required for (j> in describing an angular interval of 
2n with respect to a fixed direction is P = P^ + p( GE ) with 




p(0) 



P^ = ^ail-^+O 



(5) 
(6) 



Note that for M » m, i.e. \i -> M and v -> 0, © reduces to (JTJ (Mash- 
hoon et al 2001). Moreover, © agrees with the expression (A2.50) obtained 
by Soffel (1989) in the Brumberg representation. 



3 What Kind of Exoplanets are Suitable for Bet- 
ter Measuring the Relativistic Correction to the 
Orbital Period? 

Let us evaluate the orders of magnitude involved in a typical exoplanetary 
scenario with M = Mq, m = mj up , a = 0.1 astronomical units (AU); 
the term (37r /c 2 )^f[Ia ~ 0-1 s i while u/3 ~ 10~ 4 , so that it can safely be 
neglected in the subsequent analysis because it would be undetectable. In 
order to measure p( GE ) two requirements must be fulfilled 
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The planet must be as closest as possible to its star, i.e. a must 
be as smallest as possible. Indeed, in this case the orbital period is 
short and a large number of orbital revolutions can be recorded, thus 
increasing the observational accuracy. The observational error a P ( b B ) 
must typically be smaller than 0.1 s = 1 x 1CT 6 days. 

The Keplerian period P^ must be subtracted from the observed pe- 
riod in order to single out the relativistic component as p( GE ) = 
p(obs) _ p(o)_ As a consequence, the uncertainty in P^ due to the 
errors in the measured planet's mass and semimajor axis induces a 
systematic bias which affects the total error CTp°^ E) = o P ( bs) + cr P (o) ■ 
It turns out 



p(GE) 3 ^2 



(7) 



p(GE) 



O>(0)| a CO a 

Thus, in order to reduce the impact of a P (o) the mass of the system 
must be as largest as possible, the semimajor axis must be as smallest 
as possible and, of course, the errors in a and \i must also be small. 

In regard to the observational accuracy, the so far discovered exoplanets 
whose periods have been measured with a precision close to our stringent re- 
quirements are TrES-1 (o P ( bs) = 8x 1CP 6 days) (Alonso et al. 2004), OGLE- 
TR-56 (CTp(obB) = 5.9 x 10" 6 days) (Konacki et al. 2003) and HD209458b 
(crp(obs) = 1.8 x 10~ 7 days) (Wittenmyer et al. 2005). 



4 The HD 209458b System 

The most interesting system, for our purposes, is thus HD 209458b 'Osiris' 
which, of more than 150 recently discovered extrasolar planets, is the first 
known to transit its star (Charbonneau et al. 2000; Henry et al. 2000). 
The most recent determinations of its parameters and ephemerides, based 
on twenty-seven transit events in four bandpasses and on more than three 
years of high-precision radial velocities, are due to Wittenmyer et al. (2005). 

Its main star, 47 pc far from us, has a mass which has been estimated 
by Cody & Sasselov (2002) to lie in the range M/M© = 1.06 ± 0.13. It is 
important to note such a range does not come from a direct measurement, 
i.e. M is not one of the fifteen system's parameters simultaneously fitted 
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in the global data reduction of HD 209458 by means of the Eclipsing Light 
Curve (ELC) code (Orosz & Hauschildt 2000). It is an interval based on 
stellar evolution models and measurements of luminosity and temperature: 
no dynamical effects related in some ways to relativistic orbital motions have 
been used. M has been held fixed in (Wittenmyer et al. 2005) and different 
values for its mass have been adopted. 

The multi-color transit light curves allow to measure, among other things, 
the inclination i which amounts to 86.67 deg (Wittenmyer et al. 2005). 

The measurement of the radial velocity Km curve and of the orbital 
period P allows to determine the mass function / from which the planet's 
mass m can, thus, be deduced, for a given value of the stellar mass M. 
The planetary mass ranges from 0.590 mj up (M/M0=O.93) to 0.697 mj up 
(M/M Q =1.19) (Wittenmyer et al. 2005). Note that, up to now, the adopted 
model for f(m) is entirely Newtonian; see Section [5] for a modified version 
also including relativistic corrections. 

Finally, the semimajor axis a can be obtained from the relation which 
links it to the measured orbital period p( obs ) and the system's masses. The 
Kepler's third law © is used to this aim. The so obtained values in (Wit- 
tenmyer et al. 2005) range from 0.044 AU (M/M = 0.93) to 0.048 AU 
(M/M = 1.19). 

In regard to p( GE ), it turns out that, by using the adopted range 
for M and the measured values for a it ranges from 0.095 s to 0.112 s. 
Note that the assumption of a circular orbit for HD 209458b, adopted also 
by Wittenmyer et al. (2005) in their analysis, is consistent with a tidal 
circularization time of order 10 8 years (Bodenheimer, Laughlin & Lin 2003) 
and radial velocity observations (Mazeh et al. 2000). 

The orbital period is measured independently of any gravitational theory 
model as one of the fifteen parameters of the fit to the observations with the 
ELC software. Such important parameter has been determined with high 
accuracy by combining the ground-based observations, the transits observed 
with the photomultiplier tubes (PMTs) in the Fine Guidance Sensor (FGS) 
of the Hubble Space Telescope (HST) and the transits observed with the 
Space Telescope Imaging Spectrograph (STIS) over many cycles. The quoted 
value in (Wittenmyer et al. 2005) is 

< b i S tL my er et al.) = 3.52474554 ± 1.8 x 10- 7 days. (9) 

The uncertainty in P is thus 0.016 s. Similar results have also been obtained 
by Schultz et al. (2004) who quoted 

4°chuit z et al.) = 3.52474408 ± 2.9 x lO^days, (10) 
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with an uncertainty of 0.025 s. The situation is, thus, now more favorable 
than that described in (Iorio 2005) in which the exoplanet OGLE-TR-132b 
was examined. 

However, it must be noted that there is still a discrepancy of 0.126 s 
between the two measurements: it is as large as the relativistic effect itself. 
Such an uncertainty is due to the fact that, despite their very high pho- 
tometric precision, neither the FGS or STIS can obtain an uninterrupted 
observation of the transit due to the low orbit of the HST. The data from 
the Canadian Microvariability and Oscillations of STars (MOST) satellite, 
launched in June 2003, should overcome these problems enabling to con- 
tinuously observing HD 209458 for many transits with exceptionally high 
precision (Rucinski et al. 2003). 

The fact that the predicted relativistic component of the orbital period 
p(GE) f a \\ s j n £ m easurability domain for HD 209458b is very important 
because it opens the possibility of measuring it provided that the obtainable 
precision is high enough. As a consequence, it may turn out that relativistic 
corrections should be accounted for, e.g., in modelling both the mass func- 
tion and the orbital period in order to refine the measurements of m and a. 
In Section |S] we will deal with such topics. 



5 A Modified Mass Function to Order 0(c 



In the barycentric frame the relation 

Mcim = ma>m 



(11) 



can be assumed for the star and its planet: indeed, the post-Newtonian 
corrections of order 0(c~ 2 ) to it are negligible because are proportional to 
(Wex 1995) 



1 m 



cipn 



c 2 M(l + f) 3 



G(m + M) 



(12) 



where v and r are the relative velocity and distance. Moreover, for circular 
orbits, as in this case, (|12j) vanishes. By adding mau to both members of 
(|llj) it is possible to straightforwardly obtain 



a\[ sin 3 i 



Q 3 . '3 . 

a m sin i 



(M + m) 3 ' 
where o = ojf + o m . By using the relation 



K M 



^— p- J ™« 



in 



M + m 



2na 
~P~ 



S1I1J, 



(13) 



(14) 
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where Km is the projected semiamplitude of the star's radial velocity, is it 
possible to derive the general relation 



8tt 3 



3 '3 • ^ ■ 
a m sin i 



(m + M) 3 



(15) 



Note that the left-hand side of (|15|) is a phenomenologically determined 
quantity. 

If we use the third Kepler's law 



P 2 = 471" 



G(M + m) 



(16) 



to model the orbital period P it is possible to obtain the so called mass 
function f(M, m) as 



PK M 



m 3 sin 3 i 



2vrG (M + m) 



(17) 



For HD 209458b it amounts to 82.7 ± 1.3 m s" 1 (Wittenmyer et al. 2005). 
The relation (|17|) is of the utmost importance because it is possible to mea- 
sure the planet's mass m by keeping M fixed and measuring i, P and Km 
from the photometric transit curve and Doppler spectroscopy. 

The present-day observational accuracy in measuring the period P of 
HD 209458b suggests that general relativistic corrections might play a role 
in the system's parameters determination process. Thus, we will now derive 
a modified expression of the right-hand side of (fTTj) to order C(c~ 2 ). To this 
aim, let us write 

l3 37T 




(18) 



By defining 



4tH 



12tH 



(19) 



r 2 ; 



d 



P 1 



1)18(1 can be written to order 0(c 2 ) as 



qx 3 + bx 2 



(20) 



7 



Let us rewrite ()20j) as 



qx + b = —z; 



(21) 



since q, b, d, x are positive it admits a solution which is also unique. Let us 
look for a solution of the form 



x = x (0) (l + e), 
with e <C 1 because of order 0(c~ 2 ) and 



,(0) 



d\ 1/3 

q) 



Inserting (|22|) into (|21|) yields, to order 0{c 



3x(°)< 



Thus, 



(0) 



3<?' 



i.e. 



a ~ 



^P 2 
4^2 



2\ 1/3 



(22) 
(23) 

(24) 
(25) 
(26) 



For HD 209458b the relativistic correction to the semimajor axis amounts 
to -1 x 1(T 8 AU. 

The modified mass function can be obtained by inserting (|26|) into the 
right-hand side of (fT1)|) . By using 



a 3 ~ 



47T 2 C 2 V 27T 7 

it is straightforward to obtain 



K\ { P m 3 sin 3 i 



1 



3 /2vr/A 
c 2 V~^J 



2/3' 



2Gvr (m + M) 2 
Thus one has / = / (0) + / (GE) with, to order 0(c~ 2 
/( o) : 



(27) 



(28) 



m 3 sin 3 i 



f 



(GE) 



(M + m) 2 ' 
3m 3 sin 3 i /2vrG\ 2/3 / 1 \ 4/3 
c 2 V P ) \M + m) 



(29) 
(30) 
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For HD 209458b /(°) = 2 x 10" 10 M and /( GE ) = 
observational uncertainty in / can be evaluated as 



1.5 x l(T ib M . The 



CF f(obs) 



< 



I - — — I cr p(obs) + I 



3PK 



Al 



i x ict 11 m . 



(31) 



\2ttGJ ^ ' ' V 2vrG 

We have assumed a P ( bs) = 0.01 s (Wittenmyer et al. 2005), and <Jk m = 1-3 
m s _1 (Wittenmyer et al. 2005). This means that the relativistic correc- 
tions to the mass function are too small to play a role in the parameters' 
determination. It turns out that the major limiting factor is that due to the 
uncertainty in Km (1 x 10 -11 M ); the impact of p( obs ) is smaller amounting 
to 1 x 1O~ 17 M . 

The error in the planet's mass m, which is measured from l|17jl and is 
mainly limited by (Jk m via o - / as a m ~ Off sin 3 i, also affects the measure- 
ment of a. The uncertainty in a is 



2 

< - 
~ 3 



G(M + m) 
4vr 2 P 



nl/3 



1 

a p(obs) + - 



GP 2 



1 V3 



47r 2 (M + m) 2 



a m = 4 x 10"' AU. 

(32) 

We have used a m = 0.034 mj up (Wittenmyer et al. 2005). The impact of 
a m induces an error of 10 -7 AU. The uncertainty in the measured period 
yields an error of 10 -9 AU. Thus, the relativistic correction to a is too small 
to be measured due to o m . 



6 The Systematic Error in the Keplerian Period 

The Keplerian period must be subtracted from the measured value of 
p(° bs ) i n order to extract the relativistic correction p( GE ) , provided that the 
uncertainty in P^ is not larger than p( GE ) itself 4 . 

The precision with which m can be determined is important for our 
purposes because a m fixes, directly and indirectly via a a , the uncertainty 
in the knowledge of the Keplerian period 



a m < ^ —^—a m ^^Z^IZ aa = 4.145 s+4.157 s = 8.302 s. 

(33) 

Again, the error in m sets the limit of the measurability of the relativistic 
effect which is 75 times smaller than a P (o) . 

4 Another possible source of systematic bias is represented, in principle, by the classical 
correction P' J2 ' due to the quadrupolar mass moment J2 of the star (Iorio 2005). However, 
by assuming the fiducial value J2 ~ 10~ 6 (Miralda-Escude 2002) it turns out that p( j2 ' ~ 
-0.01 s. 
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7 Discussion and Conclusions 



In this paper we have shown that the accuracy with which it is nowadays 
possible to determine the ephemerides of the transiting extrasolar planet 
HD 209458b 'Osiris' would allow, in principle, to measure the general rel- 
ativistic correction to the orbital period P. Indeed, the orbital period of 
Osiris is measured with a ~ 0.01 s precision, while the relativistic predic- 
tion for the correction p( GE ) to the Keplerian period P(°) amounts to ~ 0.1 
s. It is likely that near-future more accurate measurements of the orbital 
periods of some of the exoplanets closest to their main stars will allow to 
extend this possibility also to them. Indeed, one of the goals of the joint 
CNES-ESA COnvection Rotation and planetary Transits (COROT) mission 
(http : //sci . esa. int/science-e/www/area/ index . cfm?f areaid=39), to 
be launched in 2006, is the discovery of new transiting giant gaseous plan- 
ets 5 close to their stars, i.e. with periods ranging from a few days to three 
months. At present, the relativistic component of the orbital periods of the 
so far discovered exoplanets, other than HD 209458b, turns out to be too 
small by one order of magnitude to be detected. 

If it was really possible to extract such a tiny slowing down of the orbital 
motion of the planet it would be the first test of general relativity, in its 
weak-field and slow-motion linearized approximation, in a planetary system 
other than our Solar System. Moreover, it would be the first measurement 
of such a consequence of the Schwarzschild gravitoelectric component of the 
space-time metric. 

Unfortunately, in the case of HD 209458b the total error <Tp°Q E) = 
o>(obs) + <7p(o) is still larger than the investigated effect. Indeed, it turns 
out that the uncertainties in the planet's mass m and semimajor axis a, due 
to the error in the projected semiamplitude of the star's radial velocity Km, 
induce a mismodelling in the Keplerian period a P (o) of 8.302 s. One-two 
orders of magnitude improvements in determining Km would be required. 
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